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ISKANDER ALIEV AND MARTIN HENK 

Abstract. Given a matrix A £ Z mxn satisfying certain regularity as- 
sumptions, a well-known integer programming problem asks to find an 
integer point in the associated knapsack polytope 

P(A, b) = {x G R> : Ax = 6} 

or determine that no such point exists. We obtain an LLL-based poly- 
nomial time algorithm that solves the problem subject to a constraint 
on the location of the vector b. 



1. Introduction and Statement of Results 

Let A G Z mxn , 1 < m < n, be an integral m x n matrix satisfying 
i) gcd (det(Ai m ) : Aj m is an m x m minor of A) = 1, 
( L1 ) ii) {x G R£ : Ax = 0} = {0}, 

where gcd(ai, . . . , a{) denotes the greatest common divisor of integers a,, 
1 < % < I. For such a matrix A and a vector b G Z m the knapsack polytope 
P(A, b) is defined as 

P(A, b) = {x € R| : ^a; = 6} . 

Observe that on account of ii), P(A, b) is indeed a polytope (or empty). 
The paper is concerned with the following integer programming problem: 

(1 2) Given input (A, b), find an integer point in P(A, b) 
or determine that no such point exists . 

The problem (j 1 . 2 [) is NP-hard (see e.g. Section 15.6 in Papadimitriou and 
Steiglitz |20j). When m = 1 we obtain the well-known integer knapsack 
problem: given integers aj, j = 1,... ,n, and b, find integers Xj > 0, j = 
1, . . . , n, such that Sj=i a j x j = b or determine that no such integers exist. 
Let us define the set 

?(A) = {b G 7L m : P(A, b)nZV 0}- 

Thus, the set F{A) will consist of all possible vectors b such that the poly- 
tope P(A, b) contains an integer point. 



2000 Mathematics Subject Classification. Primary: 90C10, 90C27, 11D07 ; Secondary: 
11H06. 

Key words and phrases. Knapsack problem; Frobenius numbers; successive minima; 
inhomogeneous minimum; distribution of lattices. 

1 



2 



ISKANDER ALIEV AND MARTIN HENK 



A set S C R m will be called a feasible set if 5 n Z m C F(A). Results of 
Aliev and Henk [2], Knight [15], Simpson and Tijdeman [26] and Pleasants, 
Ray and Simpson [21] show that the set J 7 (A) can be decomposed into 
the set of all integer points in a certain feasible (translated) cone and a 
complementary set with complex combinatorial structure. 

Note that the case m = 1 corresponds to the celebrated Frobenius problem 
and has been extensively studied in the literature. We address this problem 
below. When n = m + 1 Pleasants, Ray and Simpson [21] obtained a unique 
maximal cone whose interior is feasible. To the best of the authors knowledge 
the existence of such a maximal cone in the general case is not known. 

The location of a feasible cone is given by the diagonal Frobenius number 
defined as follows. Let Vi,...,v n € Z m be the columns of the matrix A and 
let 

C = {mv 1 H h /i„t)„ : Mi, . . . ,fJbn > 0} 

be the cone generated by v%, . . . , v n . Let also v := v\-\- . . . + v n . Following 
Aliev and Henk [2], by the diagonal Frobenius number g = g(A) of A we 
understand the minimal s > 0, such that for all b £ (sv + C) n Z m the 
polytope P(A, b) contains an integer point. Thus we have the inclusion 

(g(A)v + C)HZ m c F(A) , 

or, in other words, the translated cone g(A)v + C is feasible. 

The behavior of g{A) was investigated in Aliev and Henk [2]. The authors 
obtained an optimal up to a constant multiplier upper bound 

(1.3) g(A) < ^^-(ndet(AA T ))^ 

and estimated the expected value of the diagonal Frobenius number. 

It is natural to expect that the problem (jl.2p is solvable in polynomial 
time when the right hand side vector b belongs to a feasible cone. For such 
vectors b we a priori know that the knapsack polytope contains at least one 
integer point. We would like to propose the following conjecture. 

Conjecture 1.1. The problem is solvable in polynomial time for all 

instances (A, b) with 

b € (g(A)v + C)nZ m . 

This question is closely related to algorithmic problems in Section A.l of 
Ramirez Alfonsm |22j . 

The first result of the paper gives an estimate for the location of the 
desired feasible cone and can be considered as a step towards proving our 
conjecture. 

Theorem 1.1. There exists a polynomial time algorithm which, given (A, b), 
where A satisfies b G Z m with 

(1.4) 6 € (2( n " m )/ 2 ~ 1 p(m, n)(det( J 4,4 T )) 1 / 2 ^ + C) 
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and 

pirn, n) = 2~ 1/2 (ra - m) 1/2 n 1/2 , 
finds an integer point in the polytope P(A, b). 

In view of (|1.3p . the affirmative answer to our conjecture would imply that 
the factor 2( n ~ m ^ 2 ~ 1 p(m,n) in (jl.4p can be replaced by ("~™) TO 1 ; hence 
the exponent 2( n_m )/ 2 ~ 1 in (|1 .4[) might be redundant. 

Our next result shows that the exponent can be removed for all matrices A 
with sufficiently large det(AA T ). This phenomenon is related to the bounds 
on the efficiency of the LLL-algorithm and is a consequence of Theorem 11,41 
below. In order to state the result, let -fk be the fe-dimensional Hermite 
constant for which we refer to |18l Definition 2.2.5]. Here we just note that 
by a result of Blichfeldt (see, e.g., Gruber and Lekkerkerker 

Ik < 2 

V &k 

where at is the volume of the unit fc-ball; thus jk = 0(k). 

Theorem 1.2. There exists a polynomial time algorithm which, given (A, b), 
where A satisfies M.l\) , b E Z m with 

(1.5) b E (p(m, n)(det(AA T )) l / 2 v + C) 

and 

1.6 det(AA T ) > ± '- = !^ , 

finds an integer point in the polytope P(A, b). 

Thus, if the dimension n is concerned, Theorem 11.11 gives an exponential 
bound in n for the location of the desired feasible cone, the affirmative 
answer to Conjecture [1J] would imply the bound of order n 3 / 2 and for large 
determinants det(^4^4 T ) we obtained the bound of order n in Theorem 11.21 
In view of the size of jk, the lower bound for det(^4^4 T ) in (11. 6p has order 

n — logn+2n 

We would also like to mention an interesting consequence of Theorems 
11.11 and 11.21 The proof of Lemma 1.1 in Aliev and Henk [2] immediately 
implies that for any integer vector w in the interior int C of the cone C we 
have 

det(^ T ) \ 1/2 . 
n — m + 1/ 

It follows then from Theorem 11.11 that for every integer vector b E int C one 
can find in polynomial time an integer point in the polytope P(A, 76) for 
any integer vector 76 with 

2 (n-m)/2-l ( ) 

7 > FV ; det( AA T ) . 

n — m + 1 
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Moreover, if we assume (|1.6p to hold, then by Theorem 11.21 we can remove 
the exponential multiplier 2( ra-m )/ 2 ~ 1 from the latter inequality. 

Let us now consider the special case m = 1. Then A = a T with a = 
(osi, (Z2, • • • , a n ) T G Z n and (II. ip i) says that gcd(a) := gcd(ai, 02, • • • , a n ) = 

I. Due to the second assumption (11. ip ii) we may assume that all entries 
of a are positive. The largest integral value b such that for A = a T and 
b = (b) the polytope P(A, b) contains no integer point is called the Frobenius 
number of a , denoted by F(a). Frobenius numbers naturally appear in the 
analysis of integer programming algorithms (see, e.g., Aardal and Lenstra 
PQ, Hansen and Ryan [12], and Lee, Onn and Weismantel [E]). The general 
problem of finding F(a) has been traditionally referred to as the Frobenius 
problem. This problem is NP-hard (Ramirez Alfonsm |23} I22j ) and integer 
programming techniques are known to be an effective tool for investigating 
behavior of the Frobenius numbers, see e.g. Kannan [13], Eisenbrand and 
Shmonin [7J and Beihoffer et al [5]. 

Thus, when m = 1 the answer for the feasibility problem 

, , Given input (A, b), does the polytope P(A, b) 

^ ' ' contain an integer point? 

is affirmative for all instances (a T , b) with b > F(a). Therefore, it is natural 
to expect that for m = 1 the problem (jl.2p can be solved in polynomial time 
when b > c, for some function c = c(a). To the best of our knowledge, this 
conjecture with c = F(a) was first stated by Ramirez Alfonsm (for related 
algorithmic questions see Section A.l in |22j). Note that if the answer to 
the latter conjecture is affirmative, then the factor 2( n ~ 1 ^ 2 ~ 1 p(l, n) in (jl.4p 
can be replaced by an absolute constant. 

Let || • || denote the Euclidean norm. In the case m = 1, Theorems 

II. 11 and 11.21 deal with input instances (a T ,b), satisfying the inequalities 
b > 2( n - 1 )/ 2 - 1 p(l,n)||a||iX 1 a i and b > p(l,n)\\a\\ E" =1 <n, respectively. 
However, in this important special case, one can use slightly refined lower 
bounds for b. The bounds naturally follow from the geometric approach to 
the Frobenius problem developed in Kannan |13] and are closely related to 
the upper bound obtained in Fukshansky and Robins [8]. 

Let a[i] = (cti, . . . , Oj_i, Oj+i, • • • , a^). For m = 1 we obtain the following 
refinement of Theorems 11.11 and 11.21 

Theorem 1.3. Let 5 > 0. Then the conditions fli.^[ ) and tl.5\) in the 
statements of Theorems li.il and \1.2\ can be replaced by 

n 

(1.8) 6>2("- 1 )/ 2 - 1 (l + 5)p(l,n)^||aH||a 4 

i=l 

and 

n 

(1.9) b>(l + S)p(l,n)Y i \\a[{\\\a i , 

i=l 
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respectively. 

The proofs of Theorems 11.11 and 11.21 are based on the classical Babai's 
nearest point algorithm [4]. The algorithm is searching for a nearby lattice 
point and is built on the LLL lattice basis reduction (see Section [3]). The 
key ingredient of the proof of Theorem 11.21 is the following result. 

Theorem 1.4. Let Pk = & — 3— If L C Z n is a k-dimensional lattice 
and 61, 62, . . . ,bk is an LLL-reduced basis of L, then 

( \ 1//2 

(1.10) IIMI < (l+ r det 7U )2 ) Vndet(L), i = l,...,k. 



Note that the classical bounds for the lengths of the vectors in an LLL- 
reduced basis imply for all 1 < i < k the estimates 

INI < 2 i£ t i n 1 / 2 det(L) , 

see Lemma 13.31 below. In (jl.lOp we manage to remove the exponential mul- 
tiplier 2( fe-1 )/ 2 for integer lattices with sufficiently large determinant. 

2. Integer Knapsacks and Geometry of Numbers 

Our approach to the problem is based on Geometry of Numbers for which 
we refer to the books [BJ [Hll E] ■ 

By a lattice we will understand a discrete submodule L of a finite-dimen- 
sional Euclidean space. Here we are mainly interested in primitive lattices 
LcP, where such a lattice is called primitive if L = span K (L) n Z n . 

Recall that the Frobenius number F(a) is defined only for integer vec- 
tors a = (01, 02, • • • , a n ) with gcd(a) = 1. This is equivalent to the state- 
ment that the 1-dimensional lattice L = Za, generated by a is primitive. 
This generalizes easily to an m-dimensional lattice L C Z n generated by 
a±, . . . , a m £ Z n . Here the criterion is that L is primitive if and only if the 
greatest common divisor of all m x m-minors is 1. This is an immediate con- 
sequence of Cassels [6j Lemma 2, Chapterl] or see Schrijver [2SJ Corollary 
4.1c]. 

Hence, by our assumption (jl.ip i), the rows of the matrix A generate a 
primitive lattice La- The determinant of an m-dimensional lattice is the Tri- 
dimensional volume of the parallelepiped spanned by the vectors of a basis. 
Thus in our setting we have 

det(L A ) = ^det(AA T ). 

Now let A 6 Z mxn be a matrix satisfying the assumptions (jl.ip . By Va 
we will denote the m-dimensional subspace of W 1 spanned by the rows of A. 
The orthogonal complement of Va in I™ will be denoted as V^, so that 

Vjt = {x £ R n : Ax = 0}. 
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Furthermore, we will use the notation 

La = Va f 

for the integer sublattice contained in V^. Observe that (see e. g. [19] ) 



(2.1) det(L^) = det(L^) = ^&et{AA T ). 

For a A;-dimensional lattice L and an O-symmetric convex body K C 
span K L the ith-successive minimum of K with respect to L is defined as 

Xi(K, L) = min{A > : dim(A K n L) > i}, 1 < % < k, 

i.e., it is the smallest factor such that XK contains at least i linearly inde- 
pendent lattice points of L. 

The Minkowski's celebrated theorem on successive minima states (cf. [10, 
Theorem 23.1]) 

2& ^ 

(2.2) — det(L) < vol (K) TT UK, L) < 2 k det(L), 
k\ - LJ - 

where vol (K) denotes the volume of K. 

Let B be the unit ball in span ffi L. In the important special case K = B 
the Minkowski's theorem on successive minima can be improved (cf. |11| 
§18.4, Theorem 3]) to 

k 

(2.3) det(L) < JJ Xi(B, L) < JJ 2 det(L) . 

i=l 

3. Auxiliary results 

First, we will prove several lemmas that show a relation between the LLL 
reduction and successive minima. 

For a basis b\, 62, • • • , bfc of a lattice L in lR n we denote by b\, 62, . . . , b^ 
its Gram-Schmidt orthogonalization and by fiij the corresponding Gram- 
Schmidt coefficients, that is 

i-l 

&i = &i, b l = hi - mjbi , 2 < i < k , 
3=1 



and 



(pi, bj) 
II6-II 2 

\\ u 3\\ 



Put Xi = X(B,L), where B is the unit ball in span R L. Let us recall the 
following technical observation. 

Lemma 3.1. We have 

Xi > min \\bj\\ , i = l,2,...,k. 

j=i,i+l,...,k 
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Proof. The proof can be easily derived from the proof of Proposition 1.12 in 

cm. □ 



Recall that a lattice basis b\ , 62, ■ ■ ■ , is LLL-reduced if 

(a) |/%' I < I, for 1 < j < i < k\ 

(b) § ||6i_i|| 2 < ||6i|| 2 + ^ 2 ._ 1 ||b i _ 1 ||2 ; for 2 < i < k. 

The next lemma gives well-known upper bounds for the length the ith. 
vector of the LLL-reduced basis. 

Lemma 3.2. Suppose that the basis b%, &2> • • • , b^ is LLL-reduced. Then for 
1 < i < k the inequalities 

(3-1) ||b J || 2 <2 i - 1 ||& i || 2 

(3.2) ||6i|| 2 < 2 fc - x A 2 
hold. 

Proof. The inequalities (|3.ip and (|3.2p can be easily derived from the proofs 
of Propositions 1.6 and 1.12 in [T7j, respectively. □ 

The next result gives an upper bound for the lengths of the vectors in an 
LLL-reduced basis in terms of the determinant of the lattice. The bound is 
based on the classical estimates from Lenstra, Lenstra and Lovasz [T7] and, 
consequently, involves the exponential multiplier 2 <yk ~ 1 ^ 2 . 

Lemma 3.3. Let L C Z n be given by an LLL-reduced basis ^1,^2, . . . 
Then 

(3.3) max ||6j|| < 2 V n 1/2 det(L) . 
i=l,...,k 



Proof. By Proposition 1.12 of Lenstra, Lenstra and Lovasz [17] for any choice 
of linearly independent vectors x±, . . . , € L the inequality 

(3.4) \\bi\\ < 2V max{||xi||, . . . , ||x fc ||} 

holds. 

Put C n = [—1,1]™, i.e., C n is the n-dimensional cube of edge length 
2 centered at the origin. By a well-known result of Vaaler [27], any k- 
dimensional section of the cube C n has k- volume at least 2 k . In particular 
we have 

vol fe (C n nspan M (L)) > 2 k . 

Thus, by the Minkowski theorem on successive minima, applied to the sec- 
tion C n nspan IR (L) and L, there exist linearly independent vectors x\, . . . , x\~ G 
L such that 

1 1 £«i 1 1 00 • • • H^felloo < det(L) , 
where II • ||oo denotes the maximum norm. 
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Since X{ are nontrivial integral vectors we have 

maxjUajilloo, . . . , ||a;fc||oo} < det(L) . 
Combining the latter inequality with (j3.4j) we obtain the inequality ()3.3|) . 

□ 

The last lemma of this section gives an upper bound for the last successive 
minimum in terms of the determinant of the lattice. 

Lemma 3.4. 

Afe < \fn det(L) . 



Proof. By the Minkowski theorem on successive minima, applied to the set 
C n n span R (L) and the lattice L, and by the already mentioned result of 
Vaaler |27j . we have 

k 

JjA i (C n nspan K (L),L) < det(L) . 
i=i 

Since L C Z n , the interior of C n n span R (L) does not contain any nonzero 
point of L. This implies 

X k (C n n span R (L), L) < det(L) , 

so that 

Afc < \fn det(L) . 

□ 

4. Proof of Theorem 11.41 

If for some I < k we have ||bj|| > ||6fc|| then, similarly to the arguments 
below, it can be shown that the inequalities (|1.1U|) hold. Thus we may 
assume that is the longest vector of the basis bj., . . . , bfe. Now assume 
||bfc|| > (1 + c/ c ) 1/ ' 2 v / ndet(L) for some c k > 0. Write 



bfc = b k + ^2 Vkjbj > \^kj\ < \ ■ 



j<k 

Hence 

\\b k \\ 2 <\\b k \\ 2 + ^\\h\\ 2 - 

j<k 

Since X k < T/ndet(L) by Lemma [3~74~1 and ||bfe|| 2 < A 2 by Lemma [3TTI one 
concludes 



^ll^-ll 2 > c fc • n • (det(L)) 2 . 



4 

j<k 
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Hence there is an i < k — 1 with 



Nl 2 > 4c fc -(det(L)) 5 



For this i one obtains by (I3.2D 

Using Lemma 13.11 and (|3.1j) . one obtains \\ > \\b k \\ 2 > 2 1 " /c ||b fc || 2 > 
2 1 ~ fc (l + c k )n(det(L)) 2 > 2 1 - fe ra(det(L)) 2 . Therefore 

flA 2 >A 2 Al> C ^2 2 ( 2 - fc )(det(L)) 4 . 

Finally, Minkowski's second theorem (see (|2.3p ) implies 

< Pk 
k ~ (det(L)) 2 • 

5. The Algorithm. Proofs of Theorems 11.11 and 11.21 

5.1. Proof of Theorem 11.11 We shall now give a high level description 
of an algorithm that satisfies conditions stated in Theorem 11.11 First the 
algorithm constructs an arbitrary integer solution u to Ax = b and a ra- 
tional solution c to Ax = b with large positive coordinates. From this one 
computes an integer point z in P{A, b) as follows. One finds a close vector 
v to u — c in the lattice and considers z := u — v. The vector z is 
an integer vector, since u, v are integer. It is a solution to Ax = b, since 
u is and Av = 0. Next, observe that z = u — v = c — (c — u) — v, the 
vector (c — u) — v is short, and c has large coordinates. This will imply 
z £ P(A,b). 
Suppose that 

(5.1) b £ (/z(m, n)(det{AA T )) 1 / 2 v + C) n Z m 

with jLt(m, ri) = 2( n-m )/ 2-1 p(?7i, n). The algorithm is presented below: 

Input : (A, b) with A and 6 satisfying (jl.ip and (|5.ip respectively; 
Output : z G P(4,b)nZ"; 
Step 1 : Find a basis x\, . . . , x n _ m of and an integer solution u of the 

equation Ax = b. 
Step 2 : Find a point c G P(-A, b) with coordinates 

q > v(m,n)(det(AA T )) 1/2 , l<»<n. 

Step 3 : Apply the Babai's algorithm for finding a nearby lattice point to 
the basis £Ci, . . . , a^n-m and the point u — c. The algorithm returns 
a lattice point v £ L^. 

Step 4 : The output vector z := u — v. 
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Let us now show that the algorithm satisfies conditions of Theorem 11.11 
Step 1 can be performed in polynomial time by Corollary 5.3c of Schrijver 
[25]. 

To justify Step 2 we will need the following observation. 

Lemma 5.1. Let b 6 (tv + C) n Z m , t > 0. One can find in polynomial 
time a point c E P(A, b) with all coordinates Cj > t. 

Proof. Since b £ (tv + C), we have b = ^2 T i = i(t+5i)vi with Si > 0. Therefore 
the polytope P% = {x £ P : xi > t , 1 < 2 < ra} is not empty. By Lemma 
6.5.1 of Grotschel, Lovasz and Schrijver [9], one can find in polynomial time 
a vertex c of the polytope Pf. The point c clearly satisfies conditions of 
Lemma 15.11 □ 



On account of (|5.ip , we can apply Lemma 15.11 with a rational number 
t > fi(m,n)(det(AA T )) 1 / 2 . Thus we obtain in polynomial time a point 
c € P(A, b) with coordinates c, satisfying 

Q > n(m,n)(det(AA T )) 1/2 , l<i<n. 

The algorithm of Babai (see [1]), applied at Step 3, runs in polynomial 
time as well. Thus it is enough to show that the output vector z is in the 
polytope P(A,b). 

Clearly, the polytope P(A, b) contains a ball centered at c with radius 
r > min,- a , so that 



(5.2) r > fj,(m,n)(det(AA T )) 



1/2 



Since Au = b and v G L^, the output vector z satisfies the condition 
Az = b. Thus, by (|5.2p . it is enough to show that 

(5.3) \\z - c|| < fi(m,n)(det(AA T )) 1/2 . 
The point v, computed by Babai's algorithm, satisfies 

(5.4) ||(„_ C )_„||2 < (||6 1 || 2 + ... + ||6 n _ m || 2 )/4, 

where b±, . . . , b n _ m is an LLL-reduced basis of L\. Since \\z — c|| = — 
c) — v\\, by (|5.4p we have 

(n — m) 1 / 2 

(5.5) ||z — c||< max ||Oi||. 

2 i=l, m 

By Lemma 13.31 and the choice of fj,(m,n) we obtain the inequality (|5.3p . 



5.2. Proof of Theorem 11.21 The above algorithm satisfies the statement 
of Theorem 11.21 as well. To see this, we only need to replace /i(m, n) = 
2( n ~ m )/' 2 ~ 1 p^rn,n) by fi(m,n) = p(m,n) and to apply Theorem 11.41 instead 
of Lemma 13.31 in the end of the proof. 
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6. Case m = 1. Proof of Theorem 11.31 

First we will show that the polytope P(a T , b) contains a ball of sufficiently 
large radius whose center can be computed in polynomial time. 

Lemma 6.1. The polytope P(a T ,b) contains an (n — rn)- dimensional ball 
centered at a rational point c and of radius 

(6 1) r > b M 

The point c can be computed in polynomial time. 

Proof. The polytope P(a T , b) is the simplex with vertices = (b/ai)e{, 1 < 
i < n, where ej are the standard basis vectors. Hence the inner unit normal 
vectors of the facets of this simplex (in the hyperplane {x £ R n : a T x = 0}) 
are given by 

° ' ' e i ~ TT-779 a I ! 1 < j < «■ 



I I /III \' 3 llall 2 



Here ej denotes j-th unit vector in R n , and the facet corresponding to Uj is 
the convex hull of all vertices except (p/aA ej. 

Now let c* be the center of the maximal inscribed ball in the simplex 
P(a T ,b), and let r* be its radius. Since this maximal ball touches all facets 
of the simplex, the radius is (n — 1) times the ratio of volume to surface area. 
Standard calculations (see, e.g., Fukshansky and Robins [HI (19)], note that 
the formula contains the redundant factor l/(n — 1)) gives 



a 



E"=i ll a Hlk ' 

Furthermore, we know that for 1 < j < n, the vector c* — r* Uj has to lie in 
the facet corresponding to Uj. Hence the jth. coordinate of c* — r* uj has 
to be zero and so we find 



3 II«II2 



\\a\?) YX=x\\<mW 

The numbers c* are in general not rational. However we can find in 
polynomial time a rational approximation c of the vector c* which satisfies 
the statement of the lemma. 

□ 

Suppose that 



n 

(6.2) b> (l + ^(l,n)^||o[i]||oi, 

i=i 

where ji{m,n) = 2^ n ~ m ^ 2 ~ 1 p(m,n), as in Section 15.11 To prove Theorem 
11.31 we have to find in polynomial time an integer point in P{a T ,b). 
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Recall that a[i] = (ai, . . . , a»_i, a«+i, . . . , ajy). We propose the following 
modification of the algorithm from Section [5] for solving this problem. 
Steps 1, 3 and 4 remain the same. Step 2 will be modified as follows 

Step 2* : Find a point c such that P(a T , b) contains an (n — m)-dimensional 
ball centered at c and of radius 

b\\a\\ 

(i + *)E?=ill«Mlk' 

Step 2* is justified by Lemma [6. 11 To prove correctness of the algorithm, 
it is now enough to show that the point z obtained at Step 4 satisfies 

(6.3) \\z — c\\ < r . 

Observe that, as in the proof of Theorem II. 1[ the inequality (|5.5|) holds. 
Consequently, by Lemma 1331 and f|6. 2j) we obtain the inequality f|6.3j) . Thus 
the condition (|1.4|) can be replaced by ([1.8]) . 

Next, let us replace /i(m,n) = 2( n ~ m ^ 2 ~ 1 p(m,n) by /u(m,n) = p(m,n). 
In this case, by Theorem 11.41 (for simplicity applied with / o^./(det(L)) 2 re- 
placed by 1) we obtain the inequality (|6.3p as well. Thus the condition ([1.50 
can be replaced by ([1.90 . 
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